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Abstract. We discuss the relationship between a symmetry in the neutrino 
flavour evolution equations and neutrino flavour oscillations in the collective 
precession mode. This collective precession mode can give rise to spectral 
swaps (splits) when conditions can be approximated as homogeneous and 
isotropic. Multi-angle numerical simulations of supernova neutrino flavour 
transformation show that when this approximation breaks down, non-collective 
neutrino oscillation modes decohere kinematically, but the collective precession 
mode still is expected to stand out. We provide a criterion for significant flavour 
transformation to occur if neutrinos participate in a collective precession mode. 
This criterion can be used to understand the suppression of collective neutrino 
oscillations in anisotropic environments in the presence of a high matter density. 
This criterion is also useful in understanding the breakdown of the collective 
precession mode when neutrino densities are small. 
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1. Introduction 

Because of neutrino-neutrino forward scattering or neutrino self-interaction [1, 2, 3] 
neutrinos can experience collective flavour transformation in environments such as 
the early Universe (e.g., [4, 5, 6, 7, 8, 9]) and supernovae (e.g., [10, 11, 12, 13]) where 
neutrino number densities can be very large. This phenomenon is different from 
the conventional Mikheyev-Smirnov-Wolfenstein (MSW) effect [14, 15] in that the 
flavour evolution histories of neutrinos in collective oscillations are coupled together 
and must be solved simultaneously. The possibility and consequences of collective 
neutrino oscillations in supernovae were not well appreciated until the discovery that 
the ordinary matter can be "ignored" in such phenomena [16] and the first numerical 
demonstrations of "stepwise spectral swapping" (or "spectral split") [17, 18] which is 
the imprint left by the collective flavour transformation on neutrino energy spectra. 

Significant progress has been made towards understanding collective neutrino 
oscillations in supernovae (see, e.g., [19] for a brief review and references therein). In 
particular, Raffelt and Smirnov [20] demonstrated an adiabatic (precession) solution 
for the homogeneous, isotropic neutrino gas. This solution has been shown to agree in 
part with the results of "single-angle" simulations of supernova neutrino oscillations 
[21]. These single-angle simulations essentially neglect the anisotropic nature of the 
supernova environment by assuming that the flavour evolution histories of neutrinos 
along all trajectories are identical to those along a "representative" trajectory, usually 
taken to be the radial trajectory [10, 22]. The adiabatic precession solution requires 
that at any time all neutrinos reside in a pure collective oscillation mode, the 
"precession" mode, which, as shown by Duan et al [18, 23], would explain the spectral 
swap phenomenon in the single-angle simulations. 

However, the real supernova environment is highly inhomogencous and 
anisotropic. Here by "inhomogeneous" and "anisotropic" we are referring to the 
neutrino fields. Of course, the matter density distributions in the supernova 
environment are also likely to be inhomogeneous and anisotropic. To date there are a 
few "multi-angle" simulations [17, 18, 24, 25] which, like the single-angle calculations, 
also adopt spherically symmetric supernova models but do treat flavour evolution 
along different neutrino trajectories in a self-consistent way. These multi-angle 
calculations also exhibit spectral swaps. It is still not understood how the spectral 
swap phenomenon arises in the (anisotropic) multi-angle context. In fact, some studies 
seem to suggest that collective neutrino oscillations in the isotropic and anisotropic 
environments can be very different. For example, collective neutrino oscillations of the 
bipolar type can experience "kinematic decoherence" and be disrupted in anisotropic 
environments [26, 27]. Additionally, a very large matter background can result in 
neutrino oscillation phase differences between different neutrino trajectories in an 
anisotropic neutrino gas [10], and this effect recently has been shown to result in 
suppression of collective neutrino oscillations [28] . 

In this paper we discuss a SU(iVf) rotation symmetry in the neutrino flavour 
evolution equations, where Nf = 2 and 3 for the two-flavour and three-flavour neutrino 
mixing schemes, respectively The collective precession mode for neutrino oscillations 
can ensue from this symmetry, even in inhomogeneous, anisotropic environments. This 
result explains the puzzling observations of the spectral swapping phenomenon in both 
the single-angle and multi- angle simulations of supernova neutrino oscillations. 

The rest of this paper is organized as follows. In section 2 we lay out the 
general framework for neutrino flavour transformation and discuss the SU(iVf ) rotation 
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symmetry in the flavour evolution equations for a dense neutrino gas. In section 3 
we show how the collective precession mode for neutrino oscillations can arise from 
this symmetry in various environments. We also give criteria for when the collection 
precession mode can occur. In section 4 we present a new multi-angle simulation 
of supernova neutrino oscillations. We analyze the results of this calculation guided 
by our understanding of the collective precession mode. In section 5 we give our 
conclusions. 



2. Equations of motion and symmetries 

2.1. Neutrino flavour polarization matrix 

We are interested in collective flavour oscillations in neutrino gases in which neutrinos 
may experience only forward scattering on other particles (including other neutrinos) , 
but where no inelastic scattering occurs. When physical conditions change only slowly 
with spatial dimension, the flavour content of neutrinos can be described by semi- 
classical matrices of densities [29, 30] 

[pp(t,x)] a/3 =^2n VtP (t,x)(u a \'ip 1/tP (t,x))(ip VtP (t,x)\v0), (1) 
\Pp{t^)] a p = ^n PjP (t,x){Dp\'ipp lP (t,x)}{i} 9jP (t,x)\D ot }, (2) 

where a and f3 are flavour labels (e, (i, t), \tp v ^^ p (t,x)} is the state of a neutrino 
v (antineutrino with momentum p at time t and position 3?, ^^(p),p is the 
corresponding neutrino number density, and the summation runs over all neutrino 
(antineutrino) states. Matrices of densities defined in (1) and (2) contain two separate 
pieces of information. One is the overall number density of neutrinos or antineutrinos 
with momentum p at spacetime point (t, x): 

«p(f,«) S Wtn' (3) 

py ' \ Tr p p (t, x) ifp°<0. y ' 

Here for compactness we use 4-component vector p = [p ,p] to denote a neutrino or 
antineutrino momentum mode, where p° = \p\ for the neutrino and — \p\ for the 
antineutrino. Because neutrinos may experience only forward scattering, n p (t,x) 
satisfies the conservation equation 

(d t +p-V)n p (t,x) = 0, (4) 

where p = p/\p\ is the unit vector along the neutrino propagation direction. 

The other piece of information contained in the matrix of density is the "flavour 
polarization" of the neutrino. This is in analogy to, e.g., the spin polarization of 
an electron gas. We define a traccless "neutrino flavour polarization matrix" , or 
"polarization matrix" for short, 

P (t V \ = I n- 1 {t,x)p p (t,x)-Nf 1 \ if P °>0, 
p( ' j "l -n- 1 (t,x) Pp (t,x) + N- 1 \ if P °<0, [b) 

where I is the identity matrix in flavour space, and Nt = 2 and 3 for two-flavour 
and three-flavour mixing schemes, respectively. In (5) we define the polarization 
matrices for neutrinos and antineutrinos with opposite signs, with the understanding 
that antiparticles are "negative particles" or "holes" in the particle sea. This sign 
convention will make the equations of motion (e.o.m.) more succinct and is especially 
appropriate in the two-flavour mixing scheme where 2 and 2, the fundamental 
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representations of the SU(2) group, are equivalent [16] (see section 2.2). The e.o.m. for 
polarization matrix P p (t,x) can be derived easily from that for p p (t,x) [29, 30, 31] 
and is 

(fit + p ■ V)P p (f, x) = -i[H p (t, x), P P (t, x)]. (6) 

Throughout this paper we assume a vanishing CP-violating phase. (See [32] for a 
discussion of collective neutrino oscillations with a non-vanishing CP- violating phase.) 
The Hamiltonian for polarization matrix P p (t,x) is 

H p (t,a!) = H;f(t,x) + W v ? (rip' (t, x) , P p , (t, x)\Wp') , 

= H^ c + H matt (t,x) + H»/(n p ,(t,x), P p ,(t,x)\Vp'). (7) 

The background "neutrino field" 

H£"(ry (i, x), P p >(t, x)\V P ') = V2G F ^(1 - p ■ p')n p ,(t, x)P p ,(t, x) (8) 

p' 

is a function of both neutrino number densities n p i(t,x) and neutrino flavour 
polarization matrices P p >(t, x), and does not depend on the energy of the test neutrino. 
For convenience we will drop the symbol "Vp"' with the understanding that p 1 in 
Hp" (n p i(t 7 x) 7 P p i(t,x)) refers to all neutrino modes. In (8) we use ^ , to denote the 
integration over (p')° and p': 

/OO P 
d(p')° Up'. (9) 
p , -oo J 

The procedure implied in equation (9) is tantamount to a sum over neutrino and 
antincutrino energies and trajectory directions. The "vacuum field" 

M 2 

H ^ C = v (10) 

generates vacuum oscillations, where M is the neutrino mass matrix. Because the 
trace of a Hamiltonian has no effect on neutrino oscillations, we will take Hpt c to be 
traceless hereafter. With this convention (10) in the vacuum mass basis becomes 

where Am^ = to 2 — m 2 is the difference between the squares of the mass eigenvalues 
corresponding to mass eigenstates and \vj), and A a (a = 1, 2, . . . , 8) are the Gell- 
Mann matrices. For the two-flavour mixing scheme 

Am 2 A 3 
2p° T 

where A3 is the third Pauli matrix. The 2x2 case is analogous to the 3x3 case 
discussed above except A Q (a = 1,2,3) are the Pauli matrices. The "matter field" in 
the flavour basis in the supernova environment is 

H matt {t,x) = A(i,£c)diag[l,0,0] = V2G F n e (t, x) diag[l, 0, 0], (13) 

where Gf is the Fermi constant, and n e {t, x) is the net electron number density. 
The vacuum field \Vp c and the matter field H matt (£, x) together constitute the total 
external field, H^^i, x), which does not depend on neutrino flavours. 



HpO C = -^r^, (12) 
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2. 2. Vector representation of the polarization matrix 

An Nf x JVf , traceless, Hermitian matrix A can be written in terms of vector A as 

a 

(We have adopted the convention in [22]. We use bold-faced letters, e.g., A to denote 
3- vectors in coordinate space, sans-serif letters, e.g., A for matrices in flavour space, 
and letters with an arrow, e.g., A for vectors in flavour space.) The e.o.m. for the 
"polarization vector" P p (t,x) is [33] 

(d t +p-V)P p (t,x)=H p (t,x) x P p (t,x), (15) 

where the cross product between two vectors is defined by the structure constants f a bc 
of the SU(iV f ) group [34]: 

{Ax B) a = f abc A b B c . (16) 

The definition of the polarization vector P p given by (5) and (14) for the 
antineutrino has a different sign as compared with that in [29] and with the 8- 
dimensional Bloch vector in [33]. In addition to making the expression of H^(t,x) 
[see (8)] more compact, this convention is especially convenient in the two- flavour 
mixing scheme. In studying collective neutrino oscillations the corotating-frame 
transformation technique [16] is frequently used. A corotating-frame transformation 
corresponds to rewriting (15) in a reference frame that rotates about e*3 with angular 
frequency loo, where e^ is the unit basis vector corresponding to A3 in the vacuum 
mass basis. According to (7) and (10) this transformation is equivalent to a change 
in the momentum of the neutrino p — > p' where 

Sm 2 Sm 2 , A . „ .„_. 

W = W~ Wo p =p ' { ] 

With the traditional definition, the direction of the polarization vector P p must be 
reversed when p° changes sign under the transformation (17). The polarization vector 
defined by (5) and (14), however, is invariant under such transformations. This 
definition has already been adopted in some recent literature, e.g., [35]. 
We can define the magnitude of P p as 



iPpi = \p P \ = ] py2 < ^ { (N t - 1). (is) 

The equal sign in the above strict inequality relation applies only if the neutrino 
state is a pure (quantum) state, i.e., can be described by a single ket. In forward- 
scattering the coherence of the neutrino is not lost and, therefore, \P p (t, x)\ also obeys 
the conservation equation 

(d t +p-V)\P p (t,x)\=0. (19) 

In the two-flavour mixing scheme a notation related to the polarization vector is 
the neutrino flavour isospin (NFIS) [16]. It can be defined as 



(20) 
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if the neutrino state is a pure state. It obeys the e.o.m. 
{d t +p- V)s p (t,x) = s p (t,x) x [-H^ c - H mM {t,x) 

- 2V2G F - P ■ P') n P' (*> x )s P > (*i *)]• ( 21 ) 

p' 

Equation (21) shows that two NFIS's s p (t, x) and s p >(t, x) at the same space-time 
point are "antiferromagnetically" coupled. Note that we have defined Hamiltonian 
vector fields in the same way as [20] but differing by a minus sign from those in the 
original NFIS notation [16]. In supernovae the neutrinos in a given momentum mode 
usually are not in a pure state. Assuming that supernova neutrinos are emitted in 
pure flavour states at the neutrino sphere and subsequently encounter only forward- 
scattering, we can write 

n p (t,x)P p (t,x) = 2^ri a|P (t,x)s tt: p(t,a;), (22) 

a 

In (22) s ai p(t,x) is the NFIS that represents the flavour state of the neutrino or 
antineutrino at (i, x) which is pure v a or T> a at the neutrino sphere, and n atP (t,x) is 
the associated neutrino number density. The corresponding replacements in (21) are 

Sp ► Sp f * s a > p r , Ttp' * ^a' ,p r and ^ y ^ - 

p' a' ,p' 

Equivalently, we can insist on the definition (20) and replace n p (t,x) in (21) by 

n' p (t,x) = n p (t,x)\P(t,x)\. (23) 

In the latter approach, NFIS s p (t,x) represents the "average flavour state" of the 
neutrino, and n' p {t, x) is the "net number density" of the neutrino in the flavour state 
represented by s p (t, x). 

2.3. Symmetries and conservation laws 

From (8) it is easy to show that Hp"(n p '(t,x),P p '(t,x)) satisfies two important 
identities: 

Y,np(t 1 x)[P p (t 1 x) 1 Y\^(n p ,{t 1 x) 1 P p ,(t,x))] = Q- (24) 
p 

and 

\i{t,x)^{n pl {t,x),P pl {t,x))^{t,x) = H 1 < v (n pl (t 1 x) 1 {)(t 1 x)P pl (t 1 x)(j\t,x)), (25) 

where U(i, x) is an arbitrary unitary matrix. 
Equation (24) implies that, if 

[G,H c /(t,x)] = 0, (26) 

where G is a constant Nt x Nf traceless Hermitian matrix, then the lepton current 
L M (t, a;) with temporal and spatial components 

L°(t, x) = J2 n p (t, x) Tr [P p (t, x) G] , (27) 
p 

L(t,x) =Y,Pn P (t,x)Tr[P p (t,x)G\ (28) 
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satisfies the continuity equation 

<9 M L M (i, x) = d t L°(t, x) + V • L(t, x) = 0. 
This can be easily shown using (4), (6), (7), (24) and (26): 

M L"(f , x) = iJ2 np(t, x) Tr ([P p (t, x), H p (t, x)] G) 

p 

= i £ np(t, x) Tr ([P p (f, x), H^(i, x)] G) 
n p (i,x)Tr(P p (f,x) [H p f(t, x), G] 



(29) 



E 

p 



0. 



(30) 



Equation (25) implies that, if (26) is true, then the e.o.m. (6) for the polarization 
matrix is invariant under the global (i.e., independent of p, t and x) SU(iVf) 
transformation 



P p (t,x) 



P p (t, x) = exp(— i(/)G)Pp(t, x) exp(i^G), 



where <f> is an arbitrary constant scalar. This is because (26) implies 

exp(-i(/»G)H^ t (t,a;)exp(i ( /)G) = H^(t,x). 
Using (6), (25) and (32) it can be shown that 

(fit +p ■ V)P p (t, x) = -i te(t, x) + H v /(n p , (t, x), P p ,{t, x)), P p (t, x) 



(31) 



(32) 



(33) 



In the two-flavour mixing scheme it is obvious that A 3 (in the vacuum mass basis) 
commutes with H^o' = H^o c in (12) in the absence of ordinary matter. According to 

the above discussion, the e.o.m. (15) for all polarization vectors P p (t,x) are invariant 
under simultaneous rotations about e^. In a homogeneous and isotropic neutrino gas 
a collective neutrino oscillation mode, which is represented by the collective precession 
of all polarization vectors about e?3, can arise because of this symmetry [23, 21]. This 
precession mode of collective neutrino oscillations will ultimately cause the energy 
spectra of neutrinos with different flavours to be swapped at a critical energy E s , 
a phenomenon known as the "stepwise spectral swapping" [18]. Not surprisingly, 
the value of E s is determined by the conserved lepton number L° associated with A3 
[36, 20]. Similar conclusions have also been drawn for homogeneous, isotropic neutrino 
gases in the three-flavour mixing scheme [37, 38]. 

Approximate symmetries and conservation laws can exist for scenarios where 
matter densities are large. Noting that H matt (£,x) is invariant under any rotation 
in the v (li -v T subspace, one can diagonalize the v^-v T submatrix of H c o t (t, x) by a 



rotation (y e 



(y e , 1/^1, v T i). In this new basis the external field is written as 

x 



H p* (t, x) 



1 
1 



e/.i' 
.2 



efi' 



— diag[m.g e + 2V2p°G F n e (t, x), 



m l'r'} 



(34) 



which is approximately diagonalized for any neutrino mode p if GFn e (t,x) 3> 
|Am| J -/(2p°)| (e.g., [39]). The approximate symmetries of the neutrino system about 
A3 and As, therefore, exist in the new basis {v e , , v T i) instead of the vacuum mass 
basis [37]. 
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3. Collective precession mode for neutrino oscillations 

In this section we discuss the collective precession mode in the two-flavour mixing 
scheme. This collective mode solution can arise in various environments because of 
the symmetry discussed in section 2.3. Generalization to the full three-flavour mixing 
scheme is straightforward when the polarization matrix representation is used [37]. 



3.1. Stationary, homogeneous and isotropic environments 



First we shall use the symmetry viewpoint to discuss neutrino oscillations in the 
collective precession mode in stationary, homogeneous, isotropic environments. In 
such environments no physical quantity depends on the neutrino propagation direction 
p and neither the external field H p £* nor the neutrino number density n p o varies with 
space or time. Also in this case, the polarization matrices P p o(t) are uniform and 
isotropic, which means the neutrino self-interaction potential 

(35) 



H^(n (p0 o,P (p0 o(t)) = V2G 



F ^2 11(p')oP(p')o(t) 
(P')° 



does not depend on the momentum of the test neutrino. Suppose that the set of 
variables {f2, P p o|Vp } solve the following equations 



H' 



A 3 ~ 

'( n ( P ')°, P(p')°) + ^-^-j Pp° 



0. 



P p oA 3 



(36) 
(37) 



where ft is a scalar independent of p°, P p o are traceless Hcrmitian matrices, and L° is 
a constant. Here we adopt the the vacuum mass basis {vx,v%) if H matt (i) is negligible, 
and the flavour basis (y e , ty) (which is also the matter basis) if H matt (i) is dominant. 
Therefore, 



A3 

Hpo* ^ -(w p o - A) — 



2p° 



A 



A3 

2 



if A is negligible, 
if A is very large. 



where Am' 2 = 



jext 



[see (34)]. In cither case we have 



[H°o , A 3 ] = 0. 
Using (25), (36) and (39) we can show easily that 



P p o (t) = cxp 



i(m -</>(,) Y 



Ppo exp 



-i(nt-0o)^r 



(38) 

(39) 
(40) 



is a solution to the e.o.m. [(6) without spatial dependence]: 
d_ 



: P p o(t) 



-i [H2'+H w (n M . ) P M .(t)) 1 P p o(t)] , 



(41) 



where 4>q is a constant. The solution obtained from (36), (37) and (40) is called the 
"precession solution" . 

Equation (36) is effectively a set of 3 x N cn coupled nonlinear integral equations, 
where 



(42) 
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is the number of ncutrino/antincutrino energy modes. Raffelt and Smirnov [20] 
pointed out that (36) can be reduced to two nonlinear integral equations. This can 
be shown as follows. We define 

H p o = + H^> (p , ) o,P (p , ) o) + ^. (43) 

This Hamiltonian H p o commutes with P p o if (36) is true. This means that P p o is either 
aligned or antialigned with the vector field H p o : 

~ e„o\ P n a \ ~ 

P p o = p l p l H p o, (44) 
\H p o\ 

where e p o = +1 (—1) if P p o is aligned (antialigned) with H p a. 
Using (35), (38), and (44) one can obtain [20] 



E 



e p o n p o|P p o| 



7 y/[(u P ° A - fi)/^(nt°t) - (P 3 )] 2 + <A> : 



(45) 



V e p0 n p ou p o\P p0 \ _ = (A + Q)<ot) (46) 

'[(uy A - fi)/ M « ot ) (P 3 )] 2 + (Pi) 2 



where 



is the total neutrino number density, 

M « ot ) = V2G F nl°\ (48) 

and 

(^iEv^» (49) 

v p 

is the average polarization vector. Note that we have chosen an appropriate value of 
4>o so that (P2) = 0. Given the set of parameters {A, |P p o| = |P p o |, n p a, e p o |Vp }, (37), 
(45) and (46) can be solved for the set of quantities {fi, (Pi), (P3)} which, in turn, 

determine P p o through (44). 

It is obvious from (45) and (46) that fl has a simple dependence on A: 

0(A) = Q\ X=Q - A. (50) 

In other words, in the presence of a large matter density, | a=o can be calculated as if 
there is no ordinary matter (but with oj p o = Am /2 /2p°), and then fl can be obtained 

using (50). We note that (P) is independent of A and, therefore, in this case neutrino 
flavour transformation does not depend on the matter density except for an extra 
rotation in (40). This result is expected using the corotating- frame technique [16]. 

We note that (36) [or (44)-(46)] and (37) are not guaranteed to have a solution 
or solutions. However, if one can solve these equations, then a "precession solution" 
(40) automatically obtains because of (39). This precession solution corresponds to 
the collective precession of polarization vectors about the e*3 axis. 
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3.2. Slowly-varying, homogeneous and isotropic environments 

Once the collective precession mode discussed in section 3.1 is established in the 
homogeneous, isotropic neutrino gas, it can be expected that, as \(t) and n' ot (i) 
vary slowly with time t (but with n p o (i)/n* ot (£) fixed for all p°), the collective mode 
continues and transforms adiabatically. In other words, wc still have 



P p o (t) = cxp 



A, 



P p o (t) exp 



i0(t) 



A, 



(51) 



except that P p o (t) = P p o (n* ot (i)) has a weak dependence on t through n* ot (i). Because 
neutrinos encounter only forward-scatterings and, therefore, n p a (t) / n* ot (t) is constant, 
wc have 



at 



&(*)) 



(I 
df 



-lj^ V (()Tr(P p o(()A 3 



(52) 



(53) 



= [\(t)-<j>'(t)]n t ° t (t), (54) 



Similar to (45) and (46) we have 

6 p o n p o(t) \P p o\ 

p o ^{[^o - A(i) + V{t)]/n(t) - (h (t))} 2 + (Pi(t)) 2 

^—v € p o n p a (t ) WpO | | 

P " - aw + mvm- (mw + <aw> 2 

where /x(t) = ^(n£, ot (i)) depends on t through n^ ot (£). In (53) and (54) e p o is constant 

for the adiabatic process, and |P p o| docs not vary with time because the coherence in 
neutrino mixing is maintained. 

Using (52)-(54) (P(t)) and 4>(t) can be found. The polarization vectors P p o(t) 
can then be found by using 



where 



\H p o(t)\ 



(55) 



(56) 



The full adiabatic precession solution is then obtained using (51). 



3.3. Stationary, homogeneous but anisotropic environments 

We now consider a neutrino gas in a stationary, homogeneous environment. By 
stationary and homogeneous we mean that neither A or n p varies with space or time. 
Wc note that the environment is anisotropic if n p depends on the neutrino propagation 
direction p. We also note that P p (t, x) can vary with time and/or space even if the 
environment is stationary and homogeneous. Like in section 3.1 we assume that the 
set of quantities {tt, K, P p \ Vp} is a solution to 



Hext 
„o 



= 0, 



(57) 
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p 

^pn p Tr(p p A 3 ) = L, 
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(58) 

(59) 



where ft and L° are constant scalars, K and Z are constant vectors, and P p are 
traceless Hcrmitian matrices. Using (25), (39) and (57) we can show easily that 



P p (t, x) = exp 



i(Qt-K-x-M^ 



P p exp 



-i(SU-K-x-<h)Y 



(60) 



is a solution to the e.o.m. (6): 
(d t +p- V)P p (t,x) =i 

= -i[H^ t + H^(n^,Py(t,a ; )), P p (i,x)] , 



(Q-p-K)^-, P p (t,x) 



(61) 



where </>o is a constant. 

Like in the stationary, homogeneous and isotropic case, (57)-(59) are not 
guaranteed to have a solution or solutions. If such a solution does exist, however, the 
symmetry in the neutrino flavour evolution equations automatically gives a collective 
precession mode solution for neutrino oscillations (60). 

Equation (57) implies that P p is cither aligned or antialigncd with the vector field 



A;, 



H p ee Hp? + H^(n p ,,P p ,) + (n-p-K)^ 

= + ]T Mp-p' « ot ) (Pp) + (V-P-K 
p' 



A3 
2 : 



and 



\H 



pi 



i>\ 



where e p = ±1 for alignment and ant i- alignment, respectively, 

„tot _ 



and 



p 



(62) 

(63) 

(64) 
(65) 

(66) 



is the polarization vector averaged across the neutrino/antincutrino energy spectrum 
and depends on p. Averaging (63) over p° we obtain 



if, v 1 \ " e P n p \P p \ g. 
n - „o \HJ 



(67) 



According to (62), H p depends on (Pp), not on each individual P p . Therefore, (58), 
(59) and (67) are a closed set of (3 x N ang + 4) coupled nonlinear integral equations 
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from which we can solve for {H, K, (Pp) | Vp} given a specified set of parameters 
{A, \P p \ = \P p \,n p , e p \Wp}. Here 

^ang = £ 1 ( 68 ) 
p 

is the number of neutrino (angular) trajectories. 

As in the stationary, homogeneous and isotropic case, we are able to sum out 
the energy modes in obtaining the precession solution and, therefore, reduce the 
number of equations in the closed set by a factor of ~ N cn . Nevertheless, (67) 
can become very difficult to solve if -/V ang is more than a few. We also note that, 
in a stationary, homogeneous and anisotropic environment, for a precession solution 
(Pp(t,x)) and P p (t,x) do not necessarily lie in the same plane as they would in an 
isotropic environment. 



3.4. Slowly-varying, anisotropic environments 



Finally, we consider a neutrino gas in an environment where both X(t, x) and n p (t, x) 
vary slowly with time and/or space. Like the slowly- varying, homogeneous and 
isotropic case, we expect the collective precession mode for neutrino oscillations to 
be of the form 



P p (t, x) = exp 



-i(t>(t,x)Y 



P p (t, x) exp 



(69) 



In (69) P p (t, x) = P p (X(t, x), n v i (t, x)) has a weak dependence on time and space which 
arises from the matter density and neutrino number densities. The set of quantities 
{4>(t,x), P p (t,x))\\/p} is a solution to the following equations: 

A3 



H^o + % v (n p ,(t,x),Pp,(t,x)) - (d t +p-V)cf>(t,x)^, P p (t,x) 



0. 



d t L°(t,x) + V -L(t,x) = 0, 



where 



L°(t,x) = ^2n p (t,x)Ti P p (t,x)A 3 
p 



L(t,x) = J2pn p (t,x)Tr \p p (t,x)A 



(70) 
(71) 

(72) 

(73) 
by 



As discussed in section 3.3 equation set (70) can be reduced by a factor of N e 
summing it across the neutrino/antineutrino energy spectrum. 

For the collective precession mode we expect that all polarization vectors P p (t, x) 
precess collectively about £3. This collective precession is fully described by <fi(t,x) 

in (69), and P p (t,x) must not rotate about e^. This additional constraint makes 

(70) and (71) generally unsolvable unless all P p (t, x) lie in the same plane, just as in 
the homogeneous, isotropic case. Because (70) determines <fi(t, x) up to an arbitrary 
constant <fio, we choose an appropriate value of 4>q so that 

Tr[P p (t,x)A 2 ] =0 (74) 

for any neutrino mode p at any (t 7 x). 

It will prove to be helpful to explore static systems in more detail. In static 
systems all physical quantities including polarization vectors are independent of time 
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t. In such a system the collective precession mode (69) describes a wavy distribution 
of neutrino polarization P p (a;). If K(x) = ~V(f>(x) is constant, then P p {x) rotates 
about e*3 clockwise for a complete cycle as the neutrino travels along its world line for 
a distance of 2n/\p ■ K\. (In the normal mass hierarchy case, the polarization vector 
P p (t,x) for a neutrino with p° > rotates about e*3 counterclockwise along its world 
line.) To gain some insight into the vector K(x) we sum (70) over all neutrino modes 
and obtain 

K ( x ) ■ ^2p n p( x )P P A x ) = ^2iK x ) - v p o]n p (t,x)P pA (x). (75) 
p p 
Equation (75) shows that K(x) describes an average oscillation behaviour for 
neutrinos propagating along some average direction characteristic of the neutrino 
(lcpton) flux. The direction of K(x) can be determined easily if the system is 
fully symmetric about an axis, e.g., the z axis, at x. In this case K(x) ■ p, the 
angular precession frequency of any neutrino propagating along the direction p, must 
be the same as that for neutrinos propagating along a different direction p' as long as 
p ■ z = p' ■ z. This means that K(x) must be parallel to z. 

We note that the collective precession mode discussed here is different from the 
self-maintained coherence of neutrino oscillations in the non-spherical geometry which 
is discussed in [22]. What is proposed in [22] is based on the assumption that all 
neutrino flavour polarization vectors P p (x) are perfectly aligned or anti-aligned with 
each other. This assumption forms the basis of the single- angle approximation in the 
non-spherical geometry which allows the computation of neutrino flavour evolution 
along the "streamlines" . Here streamlines are aligned along the direction of the 
neutrino number flux 

F(x) = ^2pn p (x). (76) 

p 

In the collective precession mode, polarization vectors are not required to be aligned 
(and, in fact, cannot be perfectly aligned) with each other. The vector K{x) is 
generally not parallel to the neutrino streamlines, cither. 

Because it takes little time for neutrinos to traverse the region of high neutrino 
fluxes in supernovae, all current numerical calculations for supernova neutrino 
oscillations are carried out as if all physical conditions such as neutrino fluxes and 
the matter profile are static. Computations with various physical conditions that 
correspond to supernova evolution over time can be pieced together to give a dynamic 
picture of neutrino flavour transformation in supernovae (e.g., [40, 41]). Generally, the 
static collective mode parameters computed with physical inputs at various epochs will 
not be the same. This will give a time dependence to (j>(t, x) which actually could, at 
least in principle, be derived from (70) and (71). The dynamic collective precession 
mode (69) describes a wave-like distribution of neutrino polarization whose "phase" 
cf>(t, x) travels with velocity 

K(t,x) 
\K(t,x) 

where tt(t,x) = —<fi(t,x). It is clear that the static assumption is valid if 

\ Q \ MA m-5\ ( T dyn\~ 1 ( Am 2 \ _1 / E+ \ 



V(t,x) = n(t,x) izr ^'J , (77) 



In (78) we have taken |f2| ~ 27r/rd yn and |.RT| = Am 2 /(2£' Vr ), where Td yn is the typical 
time scale for the variation of the relevant physical conditions in supernovae. 
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3.5. Criteria for collection neutrino oscillations 

In section 3.4 we have assumed that at any space-time point the collective 
precession mode for neutrino oscillations is the same as that in a stationary, 
homogeneous environment. Therefore, the collective precession mode in a time- 
varying, inhomogeneous environment can be established only if the variation of the 
physical conditions are "gentle" , so that the collective precession mode derived from 
the stationary, homogeneous approximation in the neighborhoods of different space- 
time points can be connected smoothly. To be more specific, we require that 



Kflt + P- v)Q P (t,x)\ 

\H p {t,x)\ 

where 



« 1, (79) 



9 p (t, x) = arccos f Ej^llll ] . ( 80 ) 
V \H p (t,x)\ J 

Equation (79) is similar to the adiabatic condition used for the collective precession 
mode in homogeneous, isotropic neutrino gases [20]. Unlike the adiabatic condition 
used in the MSW mechanism, (79) (or the adiabatic condition in [20]) is applicable only 
after the adiabatic precession solution has been determined. This is because neutrinos 
themselves contribute to the total flavour evolution Hamiltonian and, therefore, help 
set the adiabatic condition. 

A more practical criterion, which leads to a necessary condition for significant 
collective neutrino oscillations to occur, can be obtained by comparing the magnitudes 
of 

H^(x) ee H^(x) - p ■ K(x)^ = [X(x) - ui p0 - p ■ K(x)}^ (81) 
and 

W/(x) = H u /(n p ,(x),P p ,(x)) = a/2G f ^(1-p-p> p <(^P P '(*)- (82) 

p' 

Here we have made the static assumption. If 

|H^(x)| » |H£»|, (83) 

then (70) requires that P p (x) is either aligned or antialigned with e*3. In other words, 
no significant flavour oscillations can occur in this case, even if neutrinos are in the 
collective precession mode. We can get a crude estimate for (81) and (82) if all the 

P p {x) are taken to be aligned or antialigned with one another and we take [see (75)] 



K{x) ~ [ YM{*)-" P °\n P {t,x)e p \P p {x)\ ^ p - P '/ '\ (84) 
p ' Y, P P n p{ x )£p\ p p{ x )\ 



in computing Hp Xt (x), where e p = +1 ( — 1) if P p (x) • e 3 > (P p (x) ■ e 3 < 0) near 
the neutrino source. If (83) is satisfied for most neutrinos in some region, then no 
significant collective neutrino oscillations arc expected to occur in that region. 

In (83), |Hp Xt (:r)| measures how big the difference is between the collective 
precession frequency along the world line of a neutrino in mode p and the intrinsic 
precession frequency of the neutrino when there is no neutrino self-interaction. On 
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the other hand, |H^(x)| measures the strength of the neutrino self- interaction which 
makes collective neutrino oscillations possible. Therefore, (83) can be intuitively 
understood as the condition under which neutrino self-interaction is not strong enough 
to entrain a neutrino mode p whose intrinsic precession frequency is too different from 
the collective one. 

The condition in (83) implies that a large matter density can suppress collective 
neutrino oscillations in anisotropic environments. This is in contrast to the expected 
behaviour in the homogeneous, isotropic case, or a supernova model where the single- 
angle approximation is employed. Matter density-driven suppression of collective 

oscillations can be understood as follows. Using (84) and assuming P p (x) to be either 
aligned or antialigned with e% (i.e., neutrinos are in pure flavour states) we obtain 

K(x) ~ V2G F n e (x)L°(x)-^^, (85) 
where we have ignored vacuum oscillation frequencies. From (82) we have 

Hp~V2G F [L (x)-p-L(x)}^. (86) 
Using (81), (85) and (86) we can rewrite (83) as 

n e (x) » |i°(a;) — p ■ L(x)\ x 

Equation (87) gives the criterion for a matter density that is large enough to suppress 
collective neutrino oscillations in the anisotropic environment. Far away from the 
neutrino source we have L°(x) ~ |£(a;)| and 

n e (x) » L°(x) ~ [ntf(x) - (x) - n^(x) + nf* (x)}, (88) 

where n*" 1 ^ <. is the number density of the neutrinos or antineutrinos that are initially 
in the a flavour state at the neutrino source. The suppression of collective neutrino 
oscillations by the large matter density in the supernova environment was first shown 
in [28]. 



1 - 



P 



^4l°(x) 



\L(x)\ 



(87) 



4. Collective neutrino oscillations in supernovae 

In section 3 we have demonstrated that the collective precession mode for neutrino 
oscillations can arise because of symmetries in the neutrino flavour evolution equations. 
There is no guarantee that the adiabatic precession solution to (70) and (71) exists 
or that the corresponding collective oscillation mode is stable. Intuitively, however, 
we do expect such a collective neutrino oscillation mode to stand out under suitable 
conditions while other non-collective oscillation modes decohere kincmatically. In 
fact, it has been shown [18] that the stepwisc-spectral-swapping phenomenon can be 
the result of collective precession of polarization vectors or neutrino flavour isospins 
(NFIS's) when the single-angle approximation is valid. In this section we present a 
new multi-angle simulation which is engineered to have collective neutrino oscillations 
occur closer to the neutrino sphere than in previous simulations. Unlike previous 
calculations, the aspects of this simulation arc more difficult to capture with a 
single-angle approximation calculation. Nevertheless, our new simulation provides 
compelling evidence for the existence of the collective precession mode close to the 
neutrino sphere. We will also highlight the qualitative, multi-angle features in the 
calculation which can be explained using the criteria discussed in section 3.5. 
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4-1- Collective precession mode 

In the new calculation we adopt the same "neutrino bulb" model as in [18] with 
the radius of the neutrino sphere taken to be R = 10 km. We also take similar 
neutrino energy spectra at the neutrino sphere, but with (E Ve ) = 10 MeV, (Ep B ) = 12 
MeV, and (E Vx ) = (Ep x ) = 13 MeV, respectively. Here \v x ) (\v x )) corresponds to 
an appropriate linear combination of and \v T ) (\v^) and \v T )). We take the 
effective vacuum mixing angle to be 6 V = 0.1, and the mass-squared difference to 
be Am 2 = —3 x 10~ 3 eV 2 (inverted neutrino mass hierarchy). Note that in the 
effective 2x2 mixing scheme which we employ, 6 V ~ 6*13 and |Am 2 | is approximately 
the atmospheric neutrino mass-squared difference Am 2 tm [42]. We adopt a simple 
analytical profile for the electron number density [13]: 



where r is the distance to the center of the proto-neutron star, and the density profile 
is parametrized by S, the entropy per baryon in units of Boltzmann's constant fee- In 
the simulation we discuss here we take S = 250. 

Because the neutrino bulb model possesses spherical symmetry, the common 
azimuthal angle (j)(r) for all NFIS's in the collective precession mode, if it exists, 
must depend only on r. As a result, the NFIS's for various neutrino trajectories must 
precess in phase about ez along r. 

In figure 1(a) we plot (si(2) (?*)}, the energy-averaged NFIS component, where 
s 1 ( 2 )(r) = s p (r) ■ e 1 ( 2 ), as a function of r for neutrinos emitted as pure v e and 
propagating along a few representative trajectories. We here do not distinguish 
between the vacuum mass basis and the flavour basis because V <C 1. Because of 
the spherical symmetry, it suffices to label different neutrino trajectories by neutrino 
emission angle which can be defined to be 



evaluated at r = R, with f being the radial direction. (Note that we use 9 and to 
denote angles in flavour space and coordinate space, respectively.) Figure 1(a) indeed 
shows that, as neutrino oscillations start at r ~ 40 km, (s(r)) begin to precess in phase 
about (=3 along r. This is a clear signature of the existence of the collective precession 
mode in our numerical result. 

In figures l(b-d) we plot si^J") for a dozen individual angular and energy bins 
for both neutrinos and antineutrinos. We observe that si(r) and S2{r) for most 
neutrinos and antineutrinos are phase-locked at low to modest radii and again at 
larger radii. The phase-lock at the low to modest radii occurs because the neutrinos 
and antineutrinos participate in the collective precession mode, and the phase-lock 
at larger radii occurs because they experience vacuum oscillations. We note that the 
precession of NFIS's in vacuum oscillations is energy dependent, and can be in the 
opposite direction to the sense of precession in the collective precession mode. This 
behaviour is especially prominent for the antineutrino modes in figure 1(d). 

Two important trends in the breakdown of the collective precession mode merit 
further elaboration. One of these trends is that along each neutrino trajectory 
antineutrinos always drop out of the collective mode earlier than neutrinos, and 
antineutrinos with smaller energies drop out earlier than those with larger energies. 




(89) 



cos -d r (p) = p ■ r 



(90) 
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Figure 1. Oscillations of sj(r) (solid lines) and S2(r) (dotted lines), the 
projection of the NFIS on e*i and e*2, respectively, as functions of r. The curves in 
the top panel are for neutrinos emitted as v e and propagating along trajectories 
with cos"#jj = 1, 0.5 and 0, respectively. Here (si(2)) represents an average over 
the initial v e energy spectrum. Clearly (s(r)) for various neutrino trajectories 
precesses about e*3 in phase along r. The curves in the bottom panels show si,2 
for individual energy and angle bins (as labelled). The neutrinos or antincutrinos 
corresponding to these curves start as pure v e or D e , as labelled. Note that 
s±(r) and S2(r) are phase-locked at low to modest radii and again at larger radii. 
This suggests that the NFIS's for both neutrinos and antineutrinos precess with 
a common frequency at low to modest radii but precess with energy-dependent 
vacuum oscillation frequencies at larger radii. In addition, at large radii, the 
NFIS's for antincutrinos precess in the opposite direction with respect to those 
for neutrinos. 



This trend can be explained using the criterion in (83). Using (75) and the spherical 
symmetry of the supernova model we obtain 

K(r) = \K(r)\ = X(r)C(r) - (w(r)), (91) 

where 

C(r) = ^(r)P P Ar) > ^ (g2) 
J2 P cos r (p) n p (r)P pA (r) 



and 



Mr)) = ^n P {r)P v Ar) 

X; p cos $ r (p) n p (r)P pA (r) 



Symmetries in collective neutrino oscillations 



18 



Combining (81) and (91) we have 

\H; x \r)\ = \X(r)[C(r) cos r - 1] - [<w(r)) cosi9 r - w p0 ]|, (94) 
~ |WpO - (w(r))| if r>£. (95) 

Note that in the inverted neutrino mass hierarchy case considered here, we have 
uj p a = Am 2 /(2p°) < for neutrinos (p° > 0) and oj p a > for antincutrinos (p° < 0). 
Because there are more neutrinos than antincutrinos, we expect (u>(r)) < in this 
calculation. From (95) one sees that |Hp Xt (r)| is generally larger for antineutrinos, and 
the smaller the energy of the antineutrino is, the larger |Hp Xt (r)| is. Therefore, along 
the same neutrino trajectory (and with the same value of |H5^( r )l)> (83) is satisfied 
at lower radii for the antincutrinos with smaller energies. These antineutrinos must 
drop out of the collective precession mode earlier than other neutrinos do. 

The second trend is that neutrinos and antincutrinos propagating along the radial 
trajectory (cos#r = 1) drop out of the collective precession mode earlier than those 
with the same energies but propagating along the tangential trajectory (cosi?_r = 0). 
This can also be explained using criterion (83). As a crude estimate of the strength 
of the neutrino self-interaction we can assume that NFIS's arc aligned or antialigncd 
with each other and, therefore, 

IHKWI <* 1 - ^1 - (f)' - \ (J)' «®<>-. (96) 

,l(5)Vco S « + ±(5) 4 ,f(f) S «l. (97) 

Clearly |H^(?*)| is the weakest along the radial trajectory (cosi9 r = cos$_r = 1) for 
which (83) is satisfied first for a given neutrino energy. 



4-2. Stepwise spectral swapping 

In figure 2 we plot survival probabilities P yv {E^ r dpi) at r = 200 km as functions of 
both neutrino energy E and emission angle for both neutrinos and antineutrinos. 
As in previous calculations, figure 2 shows that v e and v x swap their energy spectra 
at energies above ~ 8 MeV, a phenomenon termed as "stepwise spectral swapping" 
or "spectral split" . 

For the single-angle approximated supernova model or for homogeneous, isotropic 
neutrino gases the stepwise-spectral-swapping phenomenon can be explained as the 
result of the collective precession mode [18, 23, 20]. If all neutrinos in a homogeneous, 
isotropic neutrino gas remain in the collective precession mode to the very end, then 
the polarization vector becomes 

P p0 = e p o sgn(w° r - uj p o) e 3 , (98) 

where e = +1 (—1) if P p o(t) is aligned (antialigned) with H p o(t), and Wp r is the 
common angular precession velocity of P p o(t) about e*3 when — * 0. As a result, 
the neutrino survival probability becomes 

!1 for neutrinos with E < E s , 
for neutrinos with E> E s , (99) 
for antincutrinos, 
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Figure 2. The neutrino survival probabilities P UV (E, as functions of neutrino 
energy E and emission angle at r = 200 km. The left panel is for neutrinos 
and the right panel is for antincutrinos. Shading code for P vu {E,'dji) is given at 
right. 



where 



E s 



(100) 



The stepwise-spectral-swapping phenomenon in an anisotropic environment, such 
as shown in figure 2, can be understood through the collective precession mode 
discussed in section 3.4 in a way similar to above argument based on (98)-(100), 
except for the replacement — > — K\ r ^oo. In reality, however, neutrinos and 
antineutrinos are not indefinitely entrained in the collective precession mode. The 
breakdown in collective behaviour results in some interesting features as shown in 
figure 2. For example, as explained in section 4.1, antineutrinos, especially those with 
low energies and/or propagating along the radial trajectory, do not participate as well 
in the collective mode as other neutrinos do. Therefore, the simplistic prescription 
in equation (99) works better for neutrinos than for antincutrinos. In fact, according 
to figure 2, this prescription breaks down for antincutrinos with energies ~ 5 MeV 
and/or for those propagating along the radial trajectory. In addition, antineutrinos 
with energies E < 0.8 MeV do not participate in collective oscillations at all and 
are almost fully converted through the conventional matter-driven MSW mechanism. 
Antincutrinos with energies E ~ 5 MeV do participate in collective oscillations. 
However, they encounter MSW resonances after they leave the collective precession 
mode and, as a result, are (partially) converted back to their original flavours. 

Another important detail in figure 2 is that the critical energy E s is actually 
different for neutrinos propagating along different trajectories. In contrast, K (r) is 
the same along these trajectories. Two main factors contribute to the variation of 
E s from trajectory to trajectory. The first factor is that the collective precession 
frequency K(r) cos$ r along a neutrino world line depends on the neutrino propagation 
direction. It is largest along the radial trajectory and the smallest along the tangential 
trajectory. The second factor is that neutrinos propagating along different trajectories 
drop out of the collective precession mode at different radii. Neutrinos propagating 
along more radially-directed trajectories stop participating in collective oscillations 
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at smaller distances from the proto-neutron star, but neutrinos propagating along 
more tangential trajectories stay in collective oscillations out to larger distances. Our 
numerical calculations show that the collective precession frequency always decreases 
as neutrino fluxes decrease. Therefore, both contributions lead to a smaller E s for the 
radial trajectory and a larger E s for the tangential trajectory. Indeed, the left panel of 
figure 2 shows that E s for the radial and tangential trajectories differ by ~ 0.7 MeV. 

5. Conclusions 

We have demonstrated the existence of the SU(Af ) rotation symmetry in the neutrino 
flavour evolution equations. This symmetry can facilitate the establishment of the 
collective precession mode for neutrino flavour oscillations in various environments. 
The stepwise-spcctral-swapping phenomenon can develop from such a collective 
neutrino oscillation mode in, e.g., supernovae. We have also given criteria for 
significant neutrino flavour oscillations to occur if neutrinos are entrained in the 
collective precession mode. These criteria can be used to understand the suppression 
of collective neutrino oscillations in anisotropic environments in the presence of a 
large matter density [28] . These criteria also illuminate the process of the breakdown 
of collective oscillations when neutrino densities are low. The results obtained in 
both our new simulation and previous multi-angle calculations for supernova neutrino 
oscillations can be understood in terms of the collective precession mode and the 
criteria we have provided. 

There remains much to be learned about the collective precession mode for 
neutrino oscillations. This collective mode can not exist when there is no solution to 
(70) and (71). A related interesting observation is that collective oscillations, if any, 
for a symmetric system with equal numbers of neutrinos and antincutrinos are quickly 
disrupted in the presence of even an infinitesimal anisotropy and flavour equipartition 
is obtained as a result [26]. In the flavour pendulum analogy [36], the asymmetry of 
neutrinos and antineutrinos constitutes the internal spin of the flavour pendulum. The 
existence of this internal spin causes the flavour pendulum to undergo the precession 
which represents the collective precession mode for neutrino oscillations. There exists 
no collective precession mode in a symmetric neutrino system — this explains the 
finding in [26]. Meanwhile, it has been shown in [27] that a typical asymmetry in the 
neutrino and antineutrino fluxes in the supernova environment will suppress multi- 
angle decoherence and, therefore, make the collective precession mode for neutrino 
oscillations possible. 
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